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We examine phantom dark energy models produced by a field with a negative kinetic term and 
a potential that satisfies the slow roll conditions: [(l/V)(dV/d(p)] 2 < 1 and (l/V)(d 2 V/d(j> 2 ) < 1. 
Such models provide a natural mechanism to produce an equation of state parameter, w, slightly 
less than —1 at present. Using techniques previously applied to quintessence, we show that in this 
limit, all such phantom models converge to a single expression for w(a), which is a function only of 
the present-day values of £2^ and w. This expression is identical to the corresponding behavior of 
w(a) for quintessence models in the same limit. At redshifts z < 1, this limiting behavior is well fit 
by the linear parametrization, w = Wo + ra„(l — a), with w a ~ —1.5(1 + wq). 



I. INTRODUCTION 

Observations [l], indicate that approximately 70% of 
the energy density in the universe is in the form of an ex- 
otic, negative-pressure component, dubbed dark energy. 
(See Ref. Q for a recent review). Taking w to be the 
ratio of pressure to density for the dark energy: 



Pde/p 



DE, 



(1) 



recent observations suggest that w is close to — 1. For 
example, if w is assumed to be constant, then —1.1 < 
w < -0.9 [11. If we are interested in dynamical models 
for dark energy, such as those that arise from a scalar field 
<fi, then such models can be significantly simplified in the 
limit that w is close to —1. This fact was exploited in Ref. 

which examined quintessence models with a nearly 
flat potential, defined as a V(cj>) satisfying the slow-roll 
conditions: 



and 



l_dV 
V~d6 



l_fV 



« 1, 



< 1. 



(2) 



(3) 



With these assumptions, it is possible to derive a generic 
expression for w as a function of the scale factor, a, that 
provides an excellent approximation to this entire class 
of potentials. (For other approaches to this problem, see 
Refs. 0,11,0). 

Here we extend these results to the case of phantom 
models, i.e., models for which w < — 1. It was first noted 
by Caldwell 

that observational data do not rule out 
the possibility that w < — 1. These phantom dark energy 
models have several interesting properties. The density 
of the dark energy increases with increasing scale factor, 
and the phantom energy density can become infinite at a 
finite time, a condition known as the "big rip" [10l. llll . ll2l ]. 
Further, it has been suggested that the finite lifetime 
for the universe which is exhibited in these models may 



provide an explanation for the apparent coincidence be- 
tween the current values of the matter density and the 
dark energy density [IH . (See Ref. [l4[ for a related, but 
different argument). 

A simple way to achieve a phantom model is to take 
a scalar field Lagrangian with a negative kinetic term. 
Such models have well-known problems [13, EH , 

but they nonetheless provide an interesting set of rep- 
resentative phantom models, and they have been widely 
studied 0, M, HI, SI SI M, SI, Us S3, HI • Here we 
assume a negative kinetic term, and then use techniques 
similar to those in Ref. @ to derive an expression for 
w(a) for phantom models satisfying conditions and 
©. 

We assume a spatially-flat universe containing only 
nonrelativistic matter and phantom dark energy, since 
radiation can be neglected in the epoch of interest. Then 



H 1 



Pt/3, 



(4) 



where H = a/ a, a is the scale factor, px is the total en- 
ergy density, and we take h — c = 8irG — 1 throughout. 

In a phantom model with negative kinetic term and 
potential V(<p), the energy density and pressure of the 
phantom are given by 



p^~{l/2)4> 2 + V(&, 



and 



PlP = -(l/2)0 2 -V(<P), 
so that the equation of state parameter is 

(l/2)0 2 + v{4>) 
(1/2)02 _ v ^y 

The evolution of cf> is given by 

4> + 3H4> - V'(4>) = 0. 



(5) 



(6) 



(7) 



(8) 



where the prime denotes the derivative with respect to 
4>. A field evolving according to equation ([8|) rolls uphill 
in the potential. 
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This equation of motion can be rewritten as [13, 




-3(1 + w) 



1 dln(x) 
1 + - 



6 dln(a) 



(9) 



where £7^ is the density of the phantom field in units of 
the critical density (note that fi^ evolves with time). In 
equation (|9|), x = <j) 2 /2V, so that x and w are related via 



1 + w 
1 — w 



(10) 



(Note that in equation ^ we use the sign convention in 
[28| . rather than the one used in (27[). Equation @ is the 
phantom version of the quintessence equation of motion 
first derived in Ref. [29]; it differs from the quintessence 
equation only in the sign of 1 + w on the right-hand side. 

We are interested in the limit where w is near —1, so 
following Ref. [f|, we define ft = —(1+w), where we take 
ft to be small and positive. This will allow us to drop 
terms of higher order in (3. Then equation ^ becomes 



= -3/3(2 + 0) + \{2 + P)^/Wh, (11) 

where we have defined A = V'/V. The fractional density 
in phantom dark energy, fi^, evolves as 



% = 3(1 + /3)fy(l - fy). 
Combining equations (fTTj) and (JT3J) yields 



rf/? _ -3/3(2 + p) + \(2 + f3) ^/mif, 

dn^~ 3(i + /3)n^(i-n ) 



(12) 



(13) 



This can be compared to the corresponding equation for 
d*y I d£l for quintessence, where 7 = 1 + w 



d 7 -3 7 (2 - 7) + A(2 - 7)^37^0 



cifL 



3(i-7)n^(i-n ) 



(14) 



Clearly, equations (fT3j) and IT4|) predict that (3(0,,),) 
for a phantom field and 7(0^) for a quintessence field 
will evolve quite differently. Now, however, we make our 
slow-roll assumptions. As in Ref. @, we take ft -C 1 and 
we assume A to be a constant, A , given by the initial 
value of V'/V. Both of these results are a consequence 
of equations ([2]) and ([3]); see Ref. [|| for the details. Then 
equation (TT3^) becomes 



d/3 



2ft 



f2^(l — f2^ 



+ 2 A V3? 

3 "(l-n^v^" 



(15) 



This is identical to the equation one obtains, in the cor- 
responding limit, for 7(0^) for quintessence [6|. Thus, in 
this limit, l-\-w(Cl^) for quintessence and — 1 — w(Q^) for 
a phantom field evolve in exactly the same way. Using 
the exact solution for equation (fT5")) from Ref. [g], we 
obtain 




Further, when 1 + w is close to — l,we can approximate 
the solution to equation (fT2"j) as [1,0] 



0,4 = 



1 



^0 - 1 



(17) 



where fi^o is the present-day value of 0,$, and we take 
a = 1 at the present. Combining equations (HHJ) and (|17[) . 
and normalizing to w = wq at present, we obtain: 



1 + w = (l + w ) 



l)a d tanh 



^l + ^-l)^ 



1 + (O^ 1 - l)o- 




(18) 



This equation is identical to the corresponding result for 
1 + w for the case of quintessence @ • The fact that these 
two different models (phantom and quintessence) yield 
an identical form for (1 + w)/(\ + wq) is not a priori 
obvious; indeed, this result is valid only in the "slow roll" 
limit considered here. 

We now compare the approximation given in equation 
(p~8j) to exact numerical results. Kujat et al. [28j showed 
that phantom potentials with a negative kinetic term 
can be divided into three broad classes, depending on 
the late-time asymptotic behavior of V'(<j>)/V(<f)). When 
V'/V — > ±00 at late times, w^, — > —00. These models 
include, for instance, negative power laws. For models 
in which V'/V asymptotically approaches a constant, w 
also approaches a constant, w$, where wo < — 1. This 
class of models includes exponential potentials. Finally, 
when V'/V — > at late times, — + —1. This class of 
models includes, for example, positive power law poten- 
tials. This final class of models, for which w — > — 1, can 
be further subdivided into models which yield a future 
singularity (a "big rip"), and those which do not. The 
exact conditions necessary to avoid a big rip involve an 
integral function of V(4>) [28[, but for power-law poten- 
tials, the condition is simpler: for V(cf>) — cf> a , a > 4 
yiel ds a big rip, while a < 4 yields no future singularity 

mm. 

We have therefore chosen to compare our analytic ap- 
proximation for w(a) with a model from each of these 
four classes. The results are displayed in Fig. 1. As in 
the case of quintessence, the agreement between the ac- 
tual evolution and our approximation in equation (I18p is 
remarkably good, with the error in w less than 0.5%. Our 
analytic approximation is well-fit by a linear relation, (30l | 



w(a) =w + w a (l - a), 



(19) 



in the regime a > 0.5, which corresponds to a redshift 
z < 1. Further, for the value of f^o adopted here 
(f^o = 0.7), we have w a « —1.5(l+w ). This is the same 
relation between w a and wq noted empirically for linear 



measurement of the baryon acoustic oscillation (BAO) 
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FIG. 1: The evolution of w for the phantom field evolving 
in a nearly flat potential as a function of the scale factor, 
a, normalized to a = 1 at the present, with il^o = 0.7 and 
wo — —1.1. Solid curve is our analytic result for the behavior 
of phantom models with a nearly flat potential and w near 
— 1. Other curves give the true evolution for the potentials 
V(cj>) = (p 6 (dotted), V(cj>) = <t> 2 (short dash), V(cj>) = 4>~ 2 
(long dash), and V((j>) = exp(A</>) (dot-dash). 



quintessence potentials [3l|, predicted for generic slow- 
roll quintessence models [6|, and derived analytically for 
linear quintessence potentials @. (Note that the latter 
reference gives a value of —1.3 rather than —1.5, but this 
difference is not significant, as the result depends on the 
assumed value of fi^o and on the assumed definition of 
w a ; the latter is ambiguous since w[a) is not exactly lin- 
ear in a in these models). 

Given that equation (fl"8)) provides a generic predic- 
tion for the behavior of w{a) in a wide class of phantom 
models, it is useful to compare this prediction to the ob- 
servations. Further, since both quintessence and phan- 
tom models produce the same generic form for w(a) in 
the limit where V'/V is small and nearly constant, we 
can combine the results for w(a) for quintessence models 
from Ref. Q with the results for w(a) for phantom mod- 
els given here into a single likelihood plot. In Fig. 2, we 
compare this generic prediction for w(a) to the SNIa and 
baryon acoustic oscillation data. The likelihoods were 
constructed using the 60 Essence supernovae, 57 SNLS 
(Supernova Legacy Survey) and 45 nearby supernovae, 
and the data release of 30 SNe la detected by HST and 
classified as the Gold sample by Riess et al. [|| . The com- 
bined dataset can be found in Ref. m. We also add the 
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FIG. 2: The 1 — a (solid) and 2 — a (dashed) contours in 
the plane defined by the present-day values of Qj, and w (de- 
noted H^o and wo respectively) for the phantom (wo < — 1) 
or quintessence (wo > — 1) field, in the limit of a nearly flat 
potential. Blue (outer) contours are limits from SNIa alone. 
Red (inner) contours are SNIa+BAO. 



scale at zbao = 0.35 as observed by the Sloan Digital 
Sky Survey [321 ] . 

It is important to note that Fig. 2 combines two 
completely different models. The portion of the figure 
with wq > — 1 corresponds to quintessence models, while 
Wq < — 1 corresponds to phantom models. We have sim- 
ply taken advantage of the fact that both models yield 
the same form for w(a) as a function of wq and Q^q in 
the limiting case where the potential is nearly flat. 

Our main conclusion is that the analytic approxima- 
tion previous derived for quintessence with a nearly flat 
potential [6] has the same form when it is extended to 
phantom models with a negative kinetic term. Thus, 
the generic prediction for the evolution w at low redshift 
{z < 1) in such models, w(a) = wq — 1.5(1 + Wq)(1 — a), 
can be applied to both wq > —1 and wq < —1. 
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